August 28, 2012 



X 



Search for D° -> 77 at BESIII 



(N 

o 

<^ Hajime MuramatsuQ 

t> 

CN Department of Physics & Astronomy 

, , University of Rochester, Rochester, NY 14627-0171 

X 
<D 

Oh 
D 

Based on the ~2.9 fb" 1 data set taken at the nominal mass of ■0(3770), 
■^J- we report our preliminary result on a search for D° — > 77. We find no 

^ significant signal and set an upper limit on B(D° — > / y / y)/B(D° — > 7r 7r°) < 

5.8 x 10~ 3 at 90% confidence level. 
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1 Why L> u -> 77 ? 



In the Standard Model (SM), flavor-changing neutral currents (FCNC) occur only 
at the loop level, where they are highly suppressed by the GIM mechanism pQ. The 
decay of D° — > 77 must be produced by charm-changing neutral currents. From 
the short distance contributions, the decay rate for D° — > 77 is predicted to be 
3 x 10 -11 [21 El Ej. However, the long distance contributions significantly enhance 
the decay rate which is estimated to be (1 — 3) x 1CT 8 [31 Hj. This FCNC decay 
could be enhanced by new physics (NP) effects which lead to contributions at loop 
level 016]. For instance, in the framework of the Minimal Supersymmetric Standard 
Model (MSSM), the calculation shows that the decay rate for c — > wy transition could 
be 6 x 1CT 6 , which is one to two orders of magnitudes enhanced relative to the SM 
rate, by considering gluino exchange [5]. 

Experimental searches for D° — > 77 were performed by the CLEO [7] and BABAR [8] 
experiments based on data samples collected at the T(4S') peak. They found no sig- 
nificant signals. The latter experiment yields the most stringent experimental upper 
limit to date on the B(D° 77), 2.4 x 10~ 6 at 90% confidence level (C.L.). 

Here I report a preliminary result of our search for D° — > 77 based on a sample 
acquired at the BEPCII energy-symmetric e + e~ collider with the BESIII detector [9]. 
This sample was collected at a/s = 3.773 GeV with an integrated luminosity of ~2.9 
fb _1 in which the background level at this open charm threshold is expected to be 
substantially lower than that at the T(45*) peak. 

2 Analyses 

One of the major backgrounds in the analysis of D° — > 77 comes from D° — > 7r°7r°. 
I first describe our study on the decay of D° — > 7r°7r°, then move on to discuss our 
search for D° — > 77. In the end, we give our result in terms of a ratio of branching 
fractions, B(D° ^)/B(D° tt% ). 



Each 7T° candidate is constructed from a pair of good photons, each of which is defined 
as: it is based on a shower in the electromagnetic Csl crystal calorimeter (EMC) of 
the BESIII detector [9]; should not match with any charged tracks reconstructed 
by the main drift chamber; must be reconstructed within barrel (|cos#| < 0.80) or 
endcap (0.86 < |cos#| < 0.93) sections of the EMC where 9 is a polar angle with 
respect to the positron beam direction; deposited energy of the shower must be at 
least 25 (50) MeV when it is reconstructed in the barrel (endcap) section. Once a 
pair is found to be a 7r° candidate (110 < M 77 < 150 MeV/c 2 ), it is kinematically 
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constrained to the known mass [10] from which the resultant constrained momentum 
of 7T° is used for the rest of the analysis. These ir° candidates are also tested for the 
possibility that one of their showers could form a 7r° with some other shower, and are 
rejected if any other pairing were more consistent with a tt° mass. Then, we form AE 
(= E^q^o —Ebeam) for every two 7r° candidates and require (—60 < AE < 30) MeV. In 
the end, we extract our signal from a distribution of beam-constrained mass, Mb c = 
\J El eam — pd 2 c 2 where Eb ea m is the beam energy and pd is the D°(— > 7r°7r°) candidate 
momentum. 

We perform a maximum likelihood fit to the resultant Mb c distribution based on 
our data sample as shown in Figure [TJ The signal shape is represented by a double 
Gaussian and its background is described by an ARGUS background function. In 
this figure, we also overlay the expected background component based on our Monte 
Carlo (MC) samples, represented by the blue-solid histogram, which apparently fails 
to describe the observed background level in data. We attribute this discrepancy to 
the incomplete simulation of continuum process under the ^(3770) peak. 

The fit yields 4081 ± 117 signal events with 52a statistical significance. The 
significance is obtained by yj — 2 x In (C wo / ' C w ) where C wo and C w are the likelihood 
values obtained from a fit without and with the signal shapes. The \ 2 f° r the fit is 
42.0 for 56 data points (minus the 5 floated parameters). 

We estimate possible systematic uncertainties on the extracted yields that include 
fit range, background shape, the requirement on AE, and reconstruction of 7r°. The 
single largest source of systematic uncertainty is due to the assumed signal shape. 
While we use a double Gaussian to represent the signal shape, the shape of a double 
Gaussian was fixed based on our MC simulation which cannot describe the data well. 
This is currently under the investigation as of this writing. Table [T] shows a summary 
of dominant sources of the estimated systematic uncertainties. 

Table 1: Systematic uncertainties on the extracted yields of D° — > -k°7i°. 



Sources 


Rel. error (%) 


Signal shape 


±8.7 


7T° recon. 


±2.0 


MC stat. 


±0.3 


Total 


±8.9 



With the total reconstruction efficiency of (23.3 ± 0.1)%, the preliminary efficiency- 
corrected yield of D° — > 7r°7r° based on our ?/>(3770) data set is 17521 ± 500(stat.) ± 
1559(syst.) events. 
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Figure 1: A fit to the D° — > 7r°7r° candidate M& c distribution based on the ?/>(3770) 
data sample. The black points are data, the black-smooth curve represents the overall 
fit (signal plus background), and the red-dashed curve corresponds to the fitted back- 
ground shape. The blue-solid histogram represents the expected background shape 
and size based on our MC samples while the blue-dotted histogram is a fit to the data 
based on this expected MC-based background shape. 
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2.2 ZT 77 

The analysis of -D — >■ 77 starts by taking the most and the 2 nd most energetic 
photon candidates in a given event, where the photon selection criteria are the same 
as the ones for the good photons previously described except we restrict these two 
photons to be reconstructed only within the barrel section of the EMC in order to 
suppress contaminations from continuum (i.e., e + e~ — > 7* — > qq — > light hadrons, 
where q = u,d,s), including doubly radiative Bhabha events. Even though these 
photons are Doppler broadened in the lab frame, the reconstructed photon energies, 
thus including the detector resolutions, are mostly found to be at least 700 MeV. 

There are two major backgrounds in this analysis, one from the continuum pro- 
cesses as mentioned above and the other from D° — > 7r°7r°. The later poses an 
irreducible background when the two tt°s decay into four photons in which two of 
them carry away most of the initial D° momentum. Figure [2] shows the AE distri- 
bution based on our MC sample in which only generic decays of DD are present, 
and D° — > 77 is set to zero. It can be seen that the background in the signal region 
(|Ai?| < 150 MeV which corresponds to ~ ±3<r) is dominated by the events from 
D° — > 7r°7r°, represented by the blue-dotted histogram in the figure. 
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Figure 2: D° — > 77 candidate AE distribution based on MC sample in which only 
generic decays of DD are present, and D° — > 77 is set to zero. Notice that almost the 
entire background comes from two-body decays, D° — > X + tc°, denoted by the red- 
dashed histogram. In the signal region (|Ai?| < 150 MeV corresponds to ~ ±3a), the 
background is dominated by events from D° — > 7r°7r°, represented by the blue-dotted 
histogram. 
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To further suppress these backgrounds, our signal selection criteria are optimized 
based on MC samples to maximize the signal sensitivity. Each of the signal pho- 
ton candidates must have a lateral shower profile that is consistent with that of an 
isolated electromagnetic shower. We suppress photons that look like coming from 
7T° — > 77 by rejecting a candidate that forms cos 61^ > 0.6 where 6* 77 is an opening 
angle between a signal photon candidate and any other photon when the given pair 
form 110 < M 77 < 150 MeV/c 2 . Requiring each of the signal photon candidates to 
be at least 20° away from any reconstructed charged track effectively suppresses the 
doubly radiative Bhabha events. In addition, we demand all selected events to satisfy 
Eemc/p < 0.8 or Eemc/v > 1-05, where p is the momentum of the fastest recon- 
structed charged track in a given event and Eemc is the corresponding deposited 
energy in the EMC. Electrons and positons tend to give Eemc I V ~ 1. To suppress 
contaminations from the rest of the continuum processes, we require there be at least 
one charged kaon reconstructed in a given event. In the end, we extract our signals 
from a AE distribution while requiring 1860 < Mb c < 1870 MeV/c 2 which gives an 
overall reconstruction efficiency of (12.02 ± 0.1)%. 

Figure [3] shows the result of a maximum-likelihood fit to the AE distribution 
based on the ?/>(3770) data set. In the fit, the signal shape is fixed by the correspond- 
ing MC shape. The background shape consists of three parts; MC-based shape to 
represent the contamination from D° — > 7r°7r° whose size is also fixed based on our 
own observation; a 1 st order polynominal that covers the contamination from Bhabha 
events which appear smoothly over the entire AE spectrum; a I s * order exponential 
polynominal, corresponding to the rest of the backgrounds. Black points are data, the 
black-solid curve is the overall fitted curve (signal plus backgrounds), the red-dashed 
curve is the fitted total backgrounds, the green curve is a sum of the exponential and 
linear polynomials. The fit gives y 2 of 63.7 for 80 data points (minus the 4 floated 
parameters) which yields —2.9 ±7.1 signal events. This translates into an upper limit 
of 11 events at 90% confidence level (C.L.) based on the Bayesian method. 

Table [2] shows a summary of our estimation of the systematic uncertainties on the 
ratio, B(D° — > jj)/B(D° — > 7r°7r°). The uncertainties we consider include: photon 
reconstruction efficiencies; the event-wise selections that require there be at least one 
charged Kaon reconstructed and the cut on Eemc/P', the requirement on M^ c \ the 
detector resolution of the signal shape; the assumed background shape, fit range. In 
this table, sources of uncertainties, that are estimated based on the data and that are 
determined to be no more than half of the measured statistical uncertainty, are listed 
as "negligible". 

Including the estimated total systematic uncertainty, we arrive at B(D° — > jj)/B(D° — > 
7T 7r°) < 5.8 x 10~ 3 at 90% C.L. 
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Figure 3: A fit to the D° — > 77 candidate AE distribution based on the -^(3770) data 
sample. Black points are data, solid black curve is the overall fitted curve (signal 
plus backgrounds), the red-dashed curve is the fitted total background, and the green 
curve is the exponential and linear polynomials. 

Table 2: Summary of the estimated systematic uncertainties on B(D° — > 77) / B(D° — > 
7r°7r°). 



Sources 


rel. errors (%) 


Photon recon. 


5.0 


Event-wise cut 


4.9 


Cut on M bc 


2.5 


Signal resolution 


negligible 


Background shape 


negligible 


Fit range 


negligible 


MC stat. 


0.4 


Stat, from D° ir°ir° 


2.9 


Syst. from D° ir°ir° 


8.9 


Total 


12.0 
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3 Conclusions and a future prospect 

Based on the ~2.9 fb _1 data set taken at y/s = 3.773 GeV, we search for D° — > 77. We 
find no significant signal and set our preliminary upper limit on B(D° — > , yy)/B(D° — > 
n°TT°) < 5.8 x 10~ 3 at 90% C.L. With the known value of B(D° 7r°7r°) [10J, this 
corresponds to B(D° -> 77) < 4.7 x 10~ 6 . 

While we are waiting for BESIII to take more data at yfs = 3.773 GeV, there 
is an alternate analysis approach that is unique to our data sample. The produced 
■0(3770) in our sample decays into a pair of D°D°. Reconstructing one of the D° 
mesons with known exclusive modes while searching for D° — > 77 in the other D° 
decay would yield an almost background-free environment, except for the irreducible 
contamination from D° — > tt 7i° for which we have control. Such a study is also 
currently under way. 
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